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Abstract — We summarise the construction of exact axisymmetric scale- 
discretised wavelets on the sphere and on the ball. The wavelet trans- 
form on the ball relies on a novel 3D harmonic transform called the 
Fourier-Laguerre transform which combines the spherical harmonic 
transform with damped Laguerre polynomials on the radial half-line. 
The resulting wavelets, called flaglets, extract scale-dependent, spatially 
localised features in three-dimensions while treating the tangential and 
radial structures separately. Both the Fourier-Laguerre and the flaglet 
transforms are theoretically exact thanks to a novel sampling theorem on 
the ball. Our implementation of these methods is publicly available 1 1 1, 
|2| and achieves floating-point accuracy when applied to band-limited 
signals. 

I. Introduction 

Spherical wavelets have been applied successfully to numerous 
problems in astrophysics and geophysics to extract features of interest 
from signals on the sphere. But in these disciplines, signals on the 
sphere are often completed with radial information such as depth, 
redshift, or distance, in which case a full 3D analysis is required. 

II. SCALE-DISCRETISED WAVELETS ON THE SPHERE 

Scale-discretised wavelets O, (Sj allow one to probe and extract 
scale-dependent, spatially localised features in signals defined on the 
sphere. In the axisymmetric case (i.e. azimuthally symmetric when 
centered on the poles) scale-discretised wavelets reduce to needlets 
141 and are constructed through a tiling of the harmonic line, thus 
defining an exact continuous transform on the sphere. Also, both the 
spherical harmonic and the scale-discretised wavelet transforms are 
exact in the discrete setting thanks to the use of a sampling theorem 
on the sphere IS). In other words a band-limited signal, i.e. described 
by a finite number of spherical harmonic coefficients, is represented 
by a finite number of samples on the sphere without any loss of 
information. Since the wavelets are band-limited by construction 
a multiresolution algorithm is used to speed up the transform by 
capturing each wavelet scale in the minimal number of samples on the 
sphere. Our implementation of the scale-dicretised wavelet transform 
is publicly available in the S2LET package | 2 | which supports the 
C, Mat lab, IDL and Java programming languages. At present the 
code is restricted to axisymmetric wavelets but will be extended to 
directional, steerable wavelets and spin functions in a future release. 

III. Flaglets on the ball 

The starting point to construct scale-discretised wavelets on the ball 
is a novel 3D transform, the Fourier-Laguerre transform, combining 
the spherical harmonics with damped Laguerre polynomials on the 
radial half-line |1 1. We construct axisymmetric wavelets, which we 
call flaglets, by separately tiling the tangential and radial harmonic 
spaces of the Fourier-Laguerre transform. Both the Fourier-Laguerre 
and flaglet transforms are exact continuous transforms, which are 
also exact in the discrete setting thanks to a 3D sampling theorem 
on the ball. Flaglets extract scale-dependent, spatially localised an- 
gular and radial features in signals defined on the ball. Since the 
wavelets are band-limit in Fourier-Laguerre space by construction, a 
multiresolution algorithm is again introduced. Our implementations 
of these transforms on the ball are publicly available in the FLAG 
and FLAGLET packages 0]. 


IV. Applications and perspectives 

The flaglet transform probes tangential and radial structures at 
scales of interest while capturing all the information of a band- 
limited signal in the minimal number of samples on the ball. It 
is suitable for high precision analysis of 3D data that requires 
the separate treatment of angular and radial components. In future 
application we intend to exploit flaglets to study galaxy surveys, 
which are used in cosmology to study the large-scale structure of 
the Universe, specifically by confronting observations (e.g. clustering 
properties) with predictions of physical models. Galaxy surveys are 
contaminated with intrinsic uncertainties and systematics affecting the 
radial and angular dimensions differently. For example photometric 
redshifts of galaxies are estimated from colour information with 
much higher uncertainty than the estimate of the angular position 
of galaxies. Hence a separate treatment of the angular and radial 
information is needed to efficiently extract cosmological information 
from galaxy surveys and to constrain relevant physical theories. Also, 
the flaglet transform takes advantage of the sparsity of these surveys: 
gravity tends to generate a filamentary structure that is captured in a 
small number of flaglet scales, as shown in figure [T| for an N-body 
simulation. 



Fig. 1. N-body simulation (top left panel) and its flaglet coefficients for 
decreasing flaglet scales (subsequent panels from left-to-right, top-to-bottom). 
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